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Abstract 

We consider the general theoretical framework to study exclusive double diffractive 
events (EDDE). It is a powerful tool to study the picture of the pp interaction. 
Basic kinematical and dynamical properties of the process, and also normalization 
of parameters via standard processes like the exclusive vector meson production 
(EVMP), are considered in detail. As an example, calculations of the cross-sections 



in the model with three pomerons for the process p + p — >■ p + M + p are presented 
A ■ for Tevatron and LHC energies. 
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1 Itroduction 



Impressive progress of the LHC experiments stimulates new investigations in different 
areas of high energy physics. Latest weighty arguments in favor of the existence of the 
Higgs boson give us the basis for futher experiments, since we have to define exactly the 
nature of this particle. It is a problem of today to have a clear instrument to determine 
its quantum numbers and couplings to other particles. Besides this task concerning Higgs 
boson, we have to continue to study other fundamental objects of high energy physics like 
jets and particles produced in different processes. 

If we consider huge number of high energy processes which have been studied for a 
long time, we will find the one which can serve as a clear source of information about high- 
energy dynamics. It is the exclusive double diffractive event (EDDE), i.e. the process 
of the type h + h — > h* + M + h*, where h — » h* scattering is quasi-diffractive, M 
is the centrally produced particle or system of particles and "+" means large rapidity 
gap (LRG). If one takes M as a single particle produced, this is the first "genuinely" 
inelastic process which not only retains a lot of features of elastic scattering (diffractive 
patterns), but also shows clearly how the initial energy is being transformed into the 
secondary particles. General properties of such amplitudes were considered in Ref. [l]-[3]. 
Theoretical consideration of these processes on the basis of Regge theory goes back to 
papers [I]- [9]. The recent interest is related to possibly good signals of centrally produced 
Higgs bosons, heavy quarkonia, glueballs, jets, gauge bosons, system of hadrons [TO] -[29]. 

Experimental study has begun since 1970's [30] ■ As Pomerons are the driving force of 
the processes in question at high energies it is naturally to expect that glueball production 
will be favorable, if one believes that Pomerons are mostly gluonic objects. Central 
glueball production was suggested as possible origin of the total cross section rise in 
Ref. [3]. One of the early proposal for experimental investigations of centrally produced 
glueballs in EDDE was made in [31]. As to the most recent experimental studies one has 
to mention the series of results from the experiment WA102 [32]- [36]. With energy rise we 
can observe central systems with higher masses [37J, like di-jet [38J, di-gamma [39j,[40j,[41j, 
heavy quarkonia [43], di-hadron ([44] and references therein). 

EDDE gives us unique experimental possibilities for particle searches and investiga- 
tions of diffraction proper. This is due to several advantages of the process: 

• clear signature of the process: central system separated from two finally detected 
protons by LRGs (see, for example, theoretical work [45], the experimental one [46J 
and references therein); 

• possibility to use the "missing mass method" [47] that improves the mass resolution; 

• strong suppression of the background due to the J z = selection rule [IE]- [50] for 
basic processes like, for example, Higgs boson production; 

• spin-parity analysis of the central system can be done [51]- [56] to determine quantum 
numbers of the central particle; 

• interesting measurements concerning the interplay between "soft" and "hard" scales 
are possible: we can obtain basic features of the interaction region (size and shape) 
from distributions in the scattering angle (diffractive patterns) [T5] . 
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There are several proposals to realize the above properties at LHC [57]-|59j. Due to the 
complicated picture of the interaction at high luminocities (a lot of "pile-up" events) 
at the moment we only have the possibility to select EDDE without detection of final 
protons [60j,[41j. The criterium of LRGs is not sufficient for our investigations, since we 
loose basic advantages of the exclusive process. For the experiment we need special low 
luminocity runs. In fact, latest LHC experiments show that the definition of diffraction 
is rather complicated subject [HI], which needs futher investigations. 

A new physics program based on polarized proton beams and tagging of forward pro- 
tons has been launched also at STAR/RHIC [62j-[65j. These experiments can drastically 
improve our understanding of diffractive mechanisms and suppress uncertainties that can 
reach sometimes an order of magnitude. 

There are many theoretical groups now that work in this area [ID]- [29]. All these 
models need to obtain values of their parameters to make predictions for the LHC energies. 
For this purpose we can use so called "standard candle" processes, i.e. events which have 
the same theoretical ingredients for the calculations. Usually authors use the following 
processes: for high central masses (M ^> 1 GeV, perturbative mechanism of Pomeron- 
Pomeron fusion dominates) 

• 7*+p -> V+p (Exclusive Vector Meson Production, EVMP), m v > 1 GeV [66]-[68]; 

• p+p -»■ p+M+p, M = ii P]. M = 77 [59"].p].[il]. M = {QQ} (heavy quarkonia, 
Xc,b) 02], M = hh (dihadron system) [it] ; 

for low central masses (M ~ 1 GeV, nonperturbative Pomeron-Pomeron fusion) 

• 7* +p->V + p (EVMP), m v ~ 1 GeV [HHj-pl]; 

• p + p — > p + M + p, M = {qq} (light meson) or "glueball" [32]- [36], M = hh 
(dihadron system) [44] ; 

In the low-mass case the mechanism of Pomeron-Pomeron fusion is considerably nonper- 
turbative. It will be taken into account in the calculations. 

In this paper we try to consider a general model-independent framework for EDDE. 
Our model for diffractive processes is taken as an example. 

2 Exclusive vector meson photoproduction. 

In this chapter we consider EVMP, i.e. exclusive photoproduction of V = QQis states, as 
the first "standard candle" for the EDDE. Basic ingredients of the theoretical framework 
are presented in the Fig. [1] and in subsections I2.1H2.3I Below we consider perturbative 
mechanism of EVMP, i.e. the mass of a vector meson is much greater than 1 GeV in 
this case. Nonperturbative mechanisms (see Fig. [TT1) will be presented in futher works. 
Following the scheme depicted in the Fig. [T] we can write the amplitude of the process as 
a convolution of the diffractive gluon-proton amplitude T and perturbative amplitude A, 
which is considered in an appropriate approach. 
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Figure 1: Exclusive vector meson photoproduction in the NRQCD framework. W is 
the 7p collision energy, t = —A 2 is the transferred momentum squared. Directions of 
momenta are shown by arrows. Vertexes A, T and Rq are defined in the text. 

2.1 Kinematics. 

First, let us introduce the notations for four- vectors depicted in the Fig. [TJ For any 
four-vector we use two following representations, 

usual: v = (v ,v 3 ;v t ) and "light- cone" : v = {v + ,v^;v t } . (1) 

For external vectors we have: 

P= { 2W ' 2W ; °J' (2) 

f W 2 _ Q2 _ m 2 A l/2 ^2 m 2 _Q 2) X 

q= { 2W ' 2W ; °J' (3) 

= / m^ + Q^ + t W\m\ + Q 2 - 3t) + (m 2 + Q 2 ){m\ + Q 2 + t) .\ 

V 2W ' 2\ 1 / 2 (W 2 ,m 2 ,-Q 2 )W 'J' U 

where 

\(x, y, z) = x 2 + y 2 + z 2 — 2xy — 2xz — 2yz, (5) 
a = W 2 , q 2 = -Q 2 , V 2 = m 2 , A 2 = t, p 2 v = m 2 v . (6) 

Let us introduce some basic variables: 

t ^ A t ,- A ' rn 2 , Q 2 

fa = k + — , k 2 = -k + — , p =p - —q, q =q + —p, (7) 

z 

« = y (y+p' + y-cO + K ±> ( 8 ) 

A = z 2 ([1 + y Q + y A ]p' - y A q') + A ± , p v = q + A , (9) 
q' 2 = p' 2 = 0, 5 = + v 7 ^) 2 + Q 2 ™ 2 , (10) 
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™>v 4^ Q 2 2 K2 A 2 4s , . 

= y = — 2~, yg = — , t — A ± = —A , ?/a = — , yo = — , (11) 

14/ m y _ rriy m v 

4 = -- 2 = "-J + ^V-%f) ^-^(V + 4y + y_) , (12) 

s = 1 GeV 2 . 

Photon and vector meson polarization vectors in the general case (Q ^ 0) can be 
represented as follows: 



e 7± g = e 70 g = , e 7 ^_ = -e 7 2 = -1 , e 70 = ^(q + s%y Q p') , (13) 



2(vA) . . 2 
£v±Pv = £v\\Pv = , £yj_ = f j_ H (p - q ) , v £ - ■- 



ev|| = — - ^(1 - 2/a)p' + A ± ) 
my 

For high-energy photoproduction (Q — )■ 0) we have 



Q < m, m y , V-t < W, (14) 



Zm - W ZQ ~ w : Zt ~ ~W' { } 

kl ~ W5 Zv ( V+ + 2^) ' V5 ZvV ~ ] ^} ' ^ 

_*),_*,.„_;_*}, (17) 

, W . , W , — A 

f W -1 f VP - , 

p ~ p' ~ -= , 0; , g ~ g' ~ 0, -=, . (19) 



V2 ) { V2 

It is convenient to introduce two transverse polarization vectors 

e 1 = (0, 0; (1, 0)), e 2 = (0, 0; (0, 1)), (20) 
q lp, g, py, e 2 = -1, (21) 

; ; f K* COS 

e^i = -€ik 2 = €iK = -{ . I 22 
2.2 Soft and hard amplitudes. 

2.2.1 NRQCD framework for the hard part of the amplitude. 

The amplitude A of the process 7(g) + g(ki) —> V(py) + g(—k 2 ) (Fig. [2]) is calculated in 
the nonrelativistic bound state approximation (see [72]- [H] and ref. therein): 
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^ab, a/3 

b ati 



R 



r.o 



Sp \6 aP (p v - m v )e v 



rik ' 
t a f b ^_ 

ij jk V3. 



Qa/3 



(23) 
(24) 



y/3m v 

where cq is the charge of heavy quark Q, Rv$ is the absolute value of the vector meson 
radial wave function at the origin, 



4-aj-b o 



due to SU(3) group rules and 



O 



a/3 



7 



-pv-^i + fc? + i 0) {-p v k 2 + k\ + i 0) 



e 7 (fci + k 2 + ^ton, j 7 ° (k 2 + zfei^v j 7 /3 
v~Pv(fci + h) + (h + k 2 ) 2 + i 0) (-pyfc 2 + k\ + i 0) 



7 «(_fc 1 + EK±!I!v)y»(_fc 1 -fc 2 + 



2 ; 



'-Vvh + k\ + i 0) (-py(/ci + fca) + (k x + fc 2 ) 2 + i 0) 



+ (1 <-> 2). (25) 



After calculations we have 



Q/3 

£ 7 ,£y 



-lQm v \ kik 2 



Pv (4 fc 2e 7 - e^ 2 e v J +pj (e£fcie 7 - e^ey) 
+ {ey^hh + g^k 1 e v k 2 e 7 ) [k\ - k lPv + k x k 2 ) 
+ (4 e 7 fc i^2 + g^k^vhe^ {k\ - k 2 p v + k x k 2 ) 



g af} (k\ - k x p v + k x k 2 ) (kj - k 2 p v + k x k 2 ) + PyPyhh 



■ [(-Pvh + k\ + i 0) (-p v k 2 + k 2 2 +i 0) (-p v (h + k 2 ) + {h + k 2 f + i 0)] 1 . (26) 

2.2.2 Diffract ive part of the amplitude. 

Diffractive gluon-proton amplitude can be written as 

T ab, of! _ yab, cflfa = 5 ab^ (p; ^ fc^^ _ ^ fc^ (27) 
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where 



G 



a 6 



G 



a/3 



pa Tjf- 
r X r 2 



QCtfi _ gOc/3 _ 2 1 ^ = p — k 2 ± - 

-P1-P2 ' ' hfa ' ' ' hh 



(28) 



Strictly speaking, in the general case we have to write more than one transverse ten- 
sor structures, but in most cases these structures satisfy relations like the Callan-Gross 
one [75], [76]. Here we can neglect terms of the order o(t) and o(z 2 ). 
Let us denote 



e y ,ev 



OtL v daa>dpp>G a 



where 



-g aa ' + ^^(^l, a ^ 2 > a ' a '^ 2 > Q ) ~~ ^1^2, q^2, a' ~ ^1^1, a^l, 



(fcxfe) 2 - A; 2 *; 2 



(29) 



(30) 



is the tensor part of the gluon propagator in the axial gauge with the axial vecor n 
k\ + k 2 = A and 



A 



t 



n 2 0. 



We can obtain the main (diagonal) contribution for further calculations when we set 



77, 



(M) 



e y = e» (i = 1,2) 
lfc{y,y+,y-) 



{y- + z 2 m y) (y+-|(i + y) + ioj (y + + !(i + y )-io)j ',(31) 

(32) 



-1 



/ c (y, y + , y_) = y 2 - 2^ (1 + (1 + 2c 2 )y) + 4y_y + (y - Ac 2 z 2 y 2 _) + 8y 2 _y 2 + . 



In the above equation (1321) c is the corresponding trigonometric function cos0 (i = 1) 
or sin0 (i = 2). In the next subsection we will have to take the integral over <ft, which 
is equivalent to the replacement c 2 — > 1/2 that is why it is convenient to introduce the 
following function 

7T 

f(y,y+,y-) = J ^-f c (y,y+,y-) = y 2 -2z 2 v yl{l + 2y)+^y + {y-2z 2 v y 2 _) + %y 2 _y 2 + . (33) 



2.2.3 Convolution and integration. Extraction of parameters for the diffrac- 
tive amplitude. 

After all convolutions the diagonal element of the amplitude for the process 7+p — > V+p 
looks as follows 



Mi 



(2tt) 4 (kl + iO)(k 2 2 + iQ) 

7T 

'4 Txm% K, v T(i,i)T D (z v y„s) 



dy dy + dy^ 
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(27T) 



m\zl 



{yl - (y. - i 0) 2 ) 



fCyly 
1 



dy dy + dy. 



f(y,y+,y- 



T D (z v y^s) 



4 (y 2 - + zfa) (y 2 - - (y- - i 0) 2 ) (y 2 + - (y + - i 0) 2 ) ' 



— _ 16 R Vfi eQy/a^a s 



3/2 

rriy tt 



2z v 



, y+ 



(i + y) 



(34) 
(35) 

(36) 



Integral (135]) can be represented in the following form 

1 

1 



i f dy [f + X + + /_2l} , (37) 



/ T = /(£+, y-, y) = £ [2 + y 2 + r\ 1 + y\ (2 - y)] , r = sign (y+y.y) , / ± = / ±1 , (38) 



dy+ I" dy_ 

L± " / w - (*+ - 1 °) 2 ) / (y- -(y~- [ °) 2 ) 

• 0(j/±4i/ + y_), (39) 

where l/z v is replaced by infinities in the upper limits, and z m is set to zero in the 
denominator (y 2 . + z^y). These replacements do not change much the final result, but 
significantly simplify further calculations. Contribution of residues in ( 1391) is dominant, if 
we take into account only the imaginary part of the amplitude T D , and we can write in 
this case 

i 

*v*-^l -j^—dy + ^T ^). (40) 
2 J y{l + y) 2 



Now we can extract parameters for the amplitude T D of the process g + p — >■ g + p. If 
we use the model of vector dominance, the amplitude of the photoproduction looks as 

M]^ V+P = -?^K V M^ V » (41) 
oa s [m v ) 

We can fix parameters of a model for meson-proton scattering by fitting the data on 
photoproduction. Then we can use expression (|34"|) and compare it with (|4ip . Finally we 
can obtain the amplitude of gluon-proton scattering from the equality 

Xv ~ s^f) M ^ ■ (42) 

In the above expressions we do not take into account the dependence of X v on t, 
since it is rather weak. It comes from the diffractive amplitude T D = T D (s,t). It means 
that (134|) .( 135|) .( 139|) . (jlDjl and ( l4"2j) are taken at any fixed value of t. If we use the data on 
total cross-sections we have to use the following equality: 



3tt 2 



2 



dt \Ms,t)\ 2 = — / dt MT; +p ^ V+p (s,t) . (43) 
\8a s {mf,)J J 

2.3 3-Pomeron model as an example for calculations. 

In principle, we can use any model for gluon-proton scattering amplitude T . For exam- 
ple, in the case of the ordinary unintegrated gluon distribution 

T D (^) - T»(*) =► /,(*, x = ^, (44) 

2 x s 
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and the integral (1401) can be rewritten as 

s/2 



X y ~-2vr 2 I %f g (x,i?,» 2 ), (45) 



which can be finally expressed in terms of the integrated parton distribution f g (x,fi 2 ). 

2.3.1 Extraction of the model parameter in the 3-Pomeron model 

As a quantitative example, in this section we consider the case of the modified Regge- 
eikonal model with three Pomerons [77]. Model parameters with errors and fits to the 
data on J/ty photoproduction can be found in [78]. It is clear from [78] that the main 
contribution to the amplitude comes from the "Born term" of the so called "hard" ("3rd") 
Pomeron in the model. Let us write corresponding amplitudes and use (j42j) to extract 
the coupling for the gluon-proton amplitude T D : 

r^cSeW^r'", (46) 



9P \S -K 2 y 

(\ Qp 3 (t) 
i) ' (47) 



M%*» v+ > = 4l W B ° h {^) (48) 

^ 3) = ^ + ^s 2 ), (49) 

a P3 (t) = 1 + A + a% A = 0.2032 ± 0.0041, a' = 0.0937 ± 0.0029 GeV" 2 , (50) 

r pP 2 = (2.4771 ± 0.0964) GeV" 2 , r gP 2 = (2.54 ± 0.41) GeV" 2 , (51) 

= 1.11 ± 0.07 , X 2 /dof = 1.48. (52) 



From (1421) we have 



c (3) _ 3?f2 c (3) (5 o\ 



With rather good accuracy of about 1% it is possible to calculate \L 



lT ■ f dyy A (A + y) 

1 vl ~2 2 ^J (l + yKl + y/y^' 1 ' 

o 

where A = «p 3 (t) — 1. \Iy\ as a function of \t\ is depicted in the Fig. El 

Since in [78] the fitting procedure was done for the data in the interval < \t\ < 1 GeV 2 
we have to use more general equation (1431) . which leads to the replacement 



•< Iv >- 



Jd|t|e- 2B I'I \I V \ 



o 



Jrf|t|e- 2S I*I 



(55) 



\ i 
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MHA-tt'lt])] 




0.2 0.4 0.6 0.8 l.O 1 ■ 

Figure 3: Function |iy(A — a'\t\)\ versus \t\ for my = rrij/q, = 3.1 GeV (solid) and 
my = m r = 9.46 GeV (dashed). 

in ( )53l) . where B = Bq + o/ln(s/so). Finally we get the coupling value 

cfj = 6.535 ± 0.418. (56) 

Here c gp is two times larger than in [75] due to different definition and corresponds to 
3.268 ± 0.209. Further we will use the definitions of the present paper for convenience. 
Errors of c gp are estimated from the errors of all the parameters in (153]) with the replace- 
ment fl55|) . 



2.3.2 Predictions for J/ty to T production ratio. 

Coupling constant c gp have to be the same for any produced vector meson. We can use 
this fact to check the model predictions for J/^> and T mesons. For the total cross-sections 
ratio it is possible to write 

n = cx^r P (W r ) _ Ex ^- ( a s (m r ) 2 <I r >(m r ,Wr)W^ 
th ' o- 7P ->j/*p(Wj/*) ~~ r .//M.^e+ e - U,(mL) < Ij/y > (mj/tp, Wjm)Wf, 

where 

mj/q , = 3.1 GeV, T J/i5 ,^ e+e - = 5.52 ± 0.18 keV, a a (m 2 J/9 ) = 0.25, (58) 
m T = 9.46 GeV, r T _» e+e - = 1.34 ± 0.05 keV, a s (m\) = 0.182. (59) 

Results of comparison based on the experimental data [66]- [68] are presented in the 
Table [Q As you can see the model adequately describes the data. 



3 Exclusive double diffractive events. 
3.1 Kinematics. 

In the Fig. H] we illustrate the perturbative mechanism for the "bare" amplitude of the 
process p+p — > p + M+p. We use the kinematics, which corresponds to the double Regge 
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Table 1: Theoretical predictions and experimental results for the ratio 1Z at different 
values of collision energies of J/*f> and T photoproduction. Experimental data are taken 
from [66]-|68j. 



Wj /qi , GeV 


W r , GeV 




Rth. ■ 10 3 


20-30 


60-130 


4.91 ±2.23 


3.49 ± 0.64 


20-30 


130-220 


9.85 ±4.37 


4.43 ± 0.66 


20-30 


60-220 


7.21 ±2.45 


4.06 ± 1.03 


30-50 


60-130 


3.86 ± 1.55 


2.89 ±0.56 


30-50 


130-220 


7.73 ± 3.0 


3.68 ±0.59 


30-50 


60-220 


5.66 ± 1.49 


3.37 ±0.88 


50-70 


60-130 


2.87 ± 1.15 


2.47 ±0.44 


50-70 


130-220 


5.75 ±2.24 


3.13 ±0.44 


50-70 


60-220 


4.21 ± 1.12 


2.87 ±0.72 


70-90 


60-130 


2.4 ±0.99 


2.2 ±0.38 


70-90 


130-220 


4.82 ± 1.9 


2.79 ±0.37 


70-90 


60-220 


3.53 ±0.96 


2.56 ±0.63 


90-110 


60-130 


2.18 ±0.88 


2.01 ±0.34 


90-110 


130-220 


4.37 ± 1.7 


2.56 ±0.33 


90-110 


60-220 


3.2 ±0.85 


2.34 ±0.57 


110-130 


60-130 


1.85 ±0.74 


1.87 ±0.31 


110-130 


130-220 


3.7 ± 1.44 


2.38 ±0.3 


110-130 


60-220 


2.71 ±0.71 


2.18 ±0.53 


130-150 


60-130 


1.54 ±0.63 


1.76 ±0.29 


130-150 


130-220 


3.09 ± 1.23 


2.24 ±0.28 


130-150 


60-220 


2.26 ±0.63 


2.05 ±0.5 


150-170 


60-130 


1.46 ±0.61 


1.67 ±0.28 


150-170 


130-220 


2.92 ± 1.19 


2.12 ±0.27 


150-170 


60-220 


2.14 ±0.62 


1.94 ±0.47 



limit. It is convenient to use light-cone representation for momenta. The notations are 
f [s m 2 J] { m 2 [s S\ 

Pi - Pi- ^-P2 ^ I y|, 0; j , p 2 ~ p 2 - ^-Pi ^ jo, j , p\ 2 = 0, 
q = Vs + + q± = q ~' & ' 

q x = q + Ai, q 2 = -q + A 2 , A li2 ~ 6i,2Pl,2, (60) 
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Figure 4: "Bare" amplitude and kinematics of the Exclusive Double Diffractive Event 
(EDDE) in the case of perturbative Pomeron-Pomeron fusion. Absorbtion in the initial 
and final pp channels is not shown. 

where 6,2 are fractions of protons' momenta carried by gluons. From the above notations 
we can obtain the relations: 

„ A 2 2 + ^ 2 2 m 2 r _ 

t h2 = a 2 2 ~ — 1,2 ; 1,2 ~ -a 2 2 , 6,2 o 

1 — U,2 

= (gi + g 2 ) 2 ^ 66^ +t 1 + t 2 - 2^hh cos O 

Ml = 66s ~ M 2 + \t x \ + |t 2 | +2 V / M 2 "cos0 o 

, A\A 2 
COS 0o = — — - — 
|Ai||A 2 | 

s i — (Pi + <?) 2 — m2 + Q 2 + V2sq~ 

S2 = (P2 - q) 2 ~m 2 + q 2 - v / 2sg+ • (61) 

Physical region of diffractive events with two rapidity gaps is defined by the following 
kinematical cuts: 

0.01 GeV 2 < |t 1)2 | < ~ 1 GeV 2 , (62) 
M 2 

imin — —1 < 6,2 < imax ~ 0.1 , (63) 

(V^h-V^h) 2 < k <(V^h + V^h) 2 (64) 

K = 66s - M 2 < M 2 
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We can write the relations in terms of yi^ and y (rapidities of hadrons and the system M 
correspondingly). For instance: 



\y\ < yo = In 



M 

y ~ 2.5 for = 14 TeV , M = 0(100 GeV) , 

2/1,2 =77 ln ^ — t >9 65 

2 — t 1>2 

3.2 Basic ingredients. 

Let us consider basic constituents of the model framework. We begin from the case of 
perturbative representation of the Pomeron-Pomeron fusion as depicted in Fig. 2] and in 
Fig. Efe.. In the perturbative case we have to use proton-gluon amplitudes T convoluted 
with gluon-gluon fusion amplitude F (see Fig. Hj), then take into account Sudakov-like 
suppression (depicted by curved wavy lines in Fig. [5^), and, finally, absorbtive (or rescat- 
tering) corrections denoted by V-blobs. In the nonperturbative case it is not possible to 
consider Pomeron as a singlet two-gluon state, and we also have to take into account the 
interaction between the central low mass state with the protons (blobs with the dashed 
boundary in the Fig. [5b). Recently it was shown in [13] that enhanced diagrams (addi- 
tional soft interactions) can play significant role. 

3.2.1 Proton-gluon scattering. 

Diffractive part T of the amplitude is calculated the same way as in fT2"T|) : 

T «» ^ T a, {pu _ gj T au ^ T au {p ^ g> fe) ( gg } 

Here we take apart color delta-functions in all amplitudes, which give overall factor 8 
after contractions and neglect terms of the order o(£j), o(ti/m 2 ). Further we will take into 
account relations 

q a Tr = 0, ft, M lT = 0, (67) 
that allow us to replace the tensor part of gluon propagators by metric tensors. 

3.2.2 Gluon-gluon fusion. 

The main contribution to the gluon-gluon fusion tensor can be represented as 

F^u = G^Fgg^M, Gfxu — Qpuv — — — — — , (68) 

Q.iQ.2 

where F gg ^ M is related to the differential cross-section 

2 



dOgg^M _ 1 AT I 7T1 |2 1=1 



N g \Fgg^ M r = \ R ^ , (69) 



di (2 8 W + l)(2 8 f ) + l)JV? 91 " 1 ISvrM^ 
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Figure 5: Scheme of calculation of the full EDDE amplitude in the case of high (a), 
intermediate (b) and low (c) invariant mass, i.e. perturbative, intermediate and nonper- 
turbative Pomeron-Pomeron fusion correspondingly. 



where 



N 2 C -1 = 8, 2sf + 1 = 2 (q 2 = 0), \F t 



gg^M\ 



256ttM' 



da, 



gg^M 



dt 



and M = gg, QQ, 77, W + W~, ZZ or other two-particle system with final momenta 
ki t 2 and t — (qx — ki) 2 . For one particle production (M = H, G(graviton), Xb,c) Fgg^M 
can be expressed in terms of the hadronic decay width of a particle 



gg^M\ 



647T 



2 1 M-^gg 
M 



(70) 



3.2.3 Sudakov-like suppression. 

In the perturbative case there is an additional factor for the gluon-gluon fusion vertex, 
which is depicted as virtual gluon corrections in Fig. If we take into account the 
emission of virtual "soft" gluons, while prohibiting the real ones, that could fill rapidity 
gaps, this results in the Sudakov-like suppression [79J: 



exp 



dpx « s (Pt 2 



l-A 



Pt 



2tt 



J zP gg (z)dz + J^2P qg (z)dz 



(71) 
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Here 

P„(z) = 6 {1 ~ * {1 ~ * ))2 , A = I 2 = -g 2 , » = M/2, (72) 

and q is the loop momentum in Fig. HJ 

3.2.4 Soft rescattering. 

We have to take into account also the hadron-hadron "soft" interaction in the initial and 
in the final states (unitary corrections or re-scattering). It is denoted by V in Fig. [5] and 
given by the following analytical expressions: 

M u ( Pl , P2 ,A h A 2) = /ggn^,) 

x M(pi -q T ,p 2 + q T ,A 1T ,A 2T )V(s',&) , (73) 
V(s,q T ) = J rf 2 6e^yi + 2iT^ pp ( S ,6), (74) 

where A iT = Ai — — q' T , A 2 t = A 2 + qx + q' T , M. is the "bare" amplitude of the process 
p + p — > p + M + p. In the case of the eikonal representation of the elastic amplitude 



T^pp we have 



V(s,q T ) = J (fbe^e 15 ^^, (75) 



where S pp ^ pp is the eikonal function. As was shown in [56], these "outer" unitary cor- 
rections strongly reduce the value of the corresponding cross-section and change the az- 
imuthal angle dependence. The ratio 



, 2 j jd 2 A l d 2 A 2 \M u \ 



< S 2 >= ^ 4 ^ \r (76) 

jjd 2 A 1 d 2 A 2 \M\ 2 

is usually called "soft survival probability" . 

3.2.5 Full construction. Convolution and integration. EDDE luminocity. 

Let us collect all the ingredients of the framework. First of all we calculate the loop 
integral 



l q = 2tt 



d 4 q f(q + ,q- } q 2 ,. . .) 

(2k)* (? + iO) (q 2 + iO) (q 2 + iO) 



q\ + iO ~ V2s£i (q--( -4^=- ~ i0 



? 2 + i0 ~ -VSfc [ 9+ + [ -|? iO ] ), (77) 



factor 2k before the integral is introduced for convenience. It was shown in [80], that 
the leading contribution arises from the region of the integration, where momentum q is 
"Glauber-like", i.e. of the order (k + m 2 /-y/i, k_m 2 /-y/i, km), where k's are of the order 1. 
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The detailed consideration of the loop integral shows that the main contribution comes 
from the poles at qf = 



Z 



da 2 f( 2l 2l_ n 2 ) 



q " 2 4 m 2 J (#2 | 

<[t l , a [>~o v g + ^ 2 

£ 2 = g 2 + A 1 2 + 2|q1|A 1 [cos(0+^) , 

qi = g 2 + A 2 2 -2|gl|A 2 |cos(0-^). (78) 
In our case we have to replace the function / in (I78p by 

iS^(^«if^) (79) 

Finally to obtain the EDDE cross-section we can introduce the luminocity function 

1 / M 2X 2 



£-EDDE — ^9^6 \~s~ J ^ < $ >, (80) 

where for low values of A 12 we replace q\ } 2 by q in X g . Now we can write 

M 2 da EDDE A da ig-?M ff .-, 

dM z dy dQgg^M d$ gg ^ M 
For resonance production Eq. (1ST]) can be simplified to 

~^y— ~ Cbddb-^—, (82) 

where TM^ gg is the resonance width. 

As in the case of EVMP we can use any model for T D . For example, 

r°(|/R)^/ s (e?V 2 ). (83) 

Parton cross-sections and resonance widths for different gluon-gluon fusion processes 
can be found elsewhere. 

3.2.6 3-Pomeron model as an example. 

In the 3-Pomeron model the diffractive proton-gluon amplitude is similar to (|46p 

Sl -' m - qq ' ) , (84) 
So - qqi J 
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and for the EDDE luminocity we have 

4 



c 



EDDE 



(3) 

Vp 3 c 9p 1 / 8 \ 2A / 2M M 1 

2W 4sUm^ ^"Ti^^ + V 



2\ A 



\i q \ 2 < s 2 >, 



where 



(85) 



<\U\> 



- toift) (!:U-v2)^'^ fe ' F ' to '~ ) " 



-.To " 



so 



+ g -, 2/2) 2(l + A) -V^ > 4 



and the "soft survival probability" can be calculated in a simple form: 



o 

oo 

fc(r) = / bdb j Q ( bT ) e -n(s,b)-n(s>, b )- b y(8B)^ n = i(J 



16) 



where 5pp^ pp can be found in [77]. Functions < S* 2 >, |/ 9 | 2 and Cedde are presented in 

Figs ei m E 



<s z > 



0.2ft 



0.15 
0.1ft 



0.05 



2 4 6 8 10 12 14 



V7, TeV 



Figure 6: "Soft survival probability" in the case of the 3-Pomeron model for different 
values of the invariant mass: M = 30 GeV (dashed), M = 125 GeV (solid) and M = 
600 GeV (dotted). 



3.2.7 Nonperturbative mechanism of Pomeron-Pomeron fusion. Discussion. 

If the central mass produced in EDDE is low (about several GeV), it is not possible to 
use perturbative representation (see Fig. 0] and Fig. [5h) for the amplitude of the process, 
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\r„t,GeV 




Figure 7: Function \I q \ 2 versus M for different values of the parameter A in the case of 
the 3-Pomeron model. A = op 3 (0) — 1 (solid) and A = ap 3 (l) — 1 (dotted). 




Figure 8: Pomeron-Pomeron luminocity Cedde at y = (left side) and as a function of y 
for different fixed invariant masses (right side) in the case of the 3-Pomeron model. In the 
left picture solid grid line denotes M = 125 GeV and dashed one denotes M = 30 GeV. 
In the right picture curves correspond to M = 20 GeV (dashed), M = 125 GeV (solid) 
and M = 600 GeV (dotted). 

and we have to use more general "nonperturbative" form (see Fig. [5b). In this case we 
have to obtain somehow the vertex "Pomeron-Pomeron-Central system" . 

From the first principles we can write the general structure of the vertex for different 
cases [56]. For example, for the production of the system with J p (spin-parity) we have 



M 0+ 


2 ~(M| 


)2(a P (0)- 


- 1 )(/ M| + 2/ 1 ) 2 , 


(89) 


M r 


2 ~(Mi 


)2(a P (0)- 


^fohh sin 2 O , 


(90) 


M l+ 


2 ~(Mi 


)2(ap(0)- 


- 1} (JoMl + J- 1 t 1 t 2 sin 2 o + T 2 ), Fo,2 ~ o(ti), 


(91) 


M 2+ 


2 ~(M| 


)2(a P (0)- 




(92) 



with notations ( 16TT) and functions defined in [56J . The general structure of helicity am- 
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plitudes from the simple Regge behaviour was also considered in [52], [53]. Experimental 
data are in good agreement with these predictions. 

There are some attempts to obtain the vertex in special models for the Pomeron. 
In Refs. [51], [55] results were obtained from the assumption that the Pomeron acts as 
a 1 + conserved or nonconserved current. The Pomeron-Pomeron fusion based on the 
"instanton" or "glueball" dynamics was considered in [H]-[2T]. You can see also recent 
papers [22] . |23] devoted to calculations of the Pomeron-Pomeron fusion vertex in the 
nonperturbative regime. 



3.3 Standard candles for high invariant masses. 

One of the basic tasks now is to make predictions for EDDE with production of funda- 
mental particles like the Higgs boson. First of all we can check any theoretical model by 
the use of the recent data from Tevatron [37]-[l3]. "Standard candles" in the case of high 
invariant mass are exclusive central di-jet and di-gamma production. We can use also Xc,o 
production, but m Xc0 = 3.5 GeV and we have to take into account also nonperturbative 
mechanisms of the Pomeron-Pomeron fusion. 

Let us analyse the CDF data in the framework of the 3- Pomeron model. In Fig. |9]you 
can see the prediction based on the HERA data (upper dotted curve) with the value of 
Cgp presented in f )56p . As was claimed in [IT]. [12], the uncertainty on the hadronic level 
can be taken into account by rescaling of jet E?- After application of this procedure to 
our result with Exjet — 0.75Et 9 we obtain the lower dashed curve in Fig. |9l which is in 
quite a good agreement with the CDF data. 

The prediction for the exclusive central di-gamma production is shown in Fig. [JUl The 




Figure 9: CDF data on the exclusive di-jet production cross-section [38] versus the lower 
cut of Exjet and the prediction of the 3P-model. Upper dotted curve corresponds to the 
value of Cgp presented in f )56|) and obtained from fitting the HERA data on EVMP. Lower 
dashed curve is obtained by rescaling E T j et = 0.75E T)9 as was proposed in []T|,[l2]. Filled 
areas denote errors from CgJ • 
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Figure 10: CDF data on the exclusive di-gamma production cross-section [38J versus lower 

(3) 

cut of E Tn and prediction of the 3P-model. Solid curve corresponds to the value of c gp 
presented in (l56il and obtained from fitting the HERA data on EVMP. Filled areas denote 
errors from c g p . Results multiplied by 3 (dotted curve) and devided by 3 (dashed curve) 
are also shown. 

value of the predicted cross-section after CDF cuts are 

Erp > 5 GeV, \rj y \ < 1 =^ a™ d ' th = 28±8 fb, (93) 
E T > 2.5 GeV, |ry 7 | < 1 =^ a excl ' th = 0.29 ± 0.08 pfc, (94) 



which are close to the predictions [13J. As was shown in [13J, uncertainties in the gluon 
distributions can lead to a factor *jj (upper and lower curves are shown for the illustration 
of this fact in the Fig. [TUT) in the theoretical prediction. If we take into account this fact, 
CDF results [39] , SOj 



E T >5 GeV, \r) 7 \ < 1 =^ o e ^ CDF < 410 fb, (95) 
E T > 2.5 GeV, |^| < 1 =► o%* > CDF = 2.48^"^ (stat) (ay at) pb (96) 

are in agreement with these predictions. 
Our prediction for the Xc,o production is 



J^-excl, th 



dy 



= 15.9 ± 4.1 nb (97) 

y=0 



and the CDF one S3] is 

ccl, CDF 

= 76 ± 10 (stat) ± 10 (syst) nb. (98) 



rl^excl, CDF 
aa Xc,0 



dy 



Theoretical prediction is about 3-=-5 times lower. As was mentioned above, the nonpertur- 
bative Pomeron-Pomeron and Reggeon-Reggeon fusion should be taken into account since 
this is the case of an intermediate invariant mass. For example, it was shown in [82].|83j. 
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that the nonperturbative contribution can be of the same order as the perturbative one. 
For more exact estimations we have to do similar calculations in our model. 

Our approach can be checked once more in the CDF process of exclusive J/\I/ produc- 
tion p + p — > p + J/ty + p which is governed by the photon- Pomeron fusion as in EVMP. 
We can obtain the cross-section as 

= Ccdf • a 1+p ^ +p {W ) = 3.51 ± 0.45 nb, (99) 

y=0 

where Wq = \Jmj/q,^/scDF — 78 GeV and Ccdf — 5.3 • 10~ 5 (see, for example, Ref. [81] 
for details of calculations). The CDF result [13] is 

= 3.92 ± 0.25 (stat) ± 0.52 (syst) nb, (100) 

y=o 

which is in a good agreement with the prediction. The prediction ( 199]) does not depend 
much on the model (since the colliding energy W lies in the HERA interval), that is why 
it is shown only as an illustration. 

Another "standard candle" is the di-hadron production. We have not considered this 
process yet, but one can find some results in |44j,[84j. 



da 



excl, th 
J/*I> 



dy 



i excL CDF 



dy 



3.4 Standard candles for low invariant masses. 

It is shown in the previous section that for invariant masses about several GeV (interme- 
diate case) our model gives results systematically lower than the experimental data. It is 
naturally to assume that we have to take into account nonperturbative contributions to 
the Pomeron-Pomeron fusion. We can use, for example, NRQCD as in the case of EVMP, 
or the method considered in [83] . where authors represent the Pomeron-Pomeron fusion 
in a similar way as the gamma-gamma fusion with different coupling (see Fig. [5b). We 
could use a similar representation in the EVMP process as depicted in Fig. [TTh to obtain 
parameters of the model. 




a) b) 

Figure 11: Representations for EVMP processes with intermediate (a) M ~ 3 — 10 GeV) 
and low (b) M ~ 1 GeV) mass of vector meson. 
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The case of M ~ 1 GeV is more complicated, because we can not use any (even semi- 
perturbative like NRQCD) mechanism for calculations. We can only restrict ourselves by 
general representations for amplitudes like (j59|) - (P?2]) and fit the data [22]-[3E] with these 
distributions. The production of low mass resonances was considered in [56]. We can 
use this "standard candle" and also EVMP process with light vector mesons depicted in 
the Fig. [TTb to obtain parameters of any nonperturbative approach. This case is more 
convenient from the experimental point of view, since cross-sections are much larger than 
for high invariant masses. This is a powerful tool to look for new states like "glueballs" 
by the use of the azimuthal angular distributions [5l].[55]. 

3.5 LHC predictions in the 3-Pomeron model. 

In this section we collect some predictions of the 3-Pomeron model for the LHC. 
The first prediction is devoted to the SM Higgs boson production 

(t p+ p-+ p+ h+p{M h = 125 GeV) ~ 0.55 ± 0.15 f&, (101) 
10" 4 < 6,2 < 0.1, 0.001 GeV 2 < |ti, 2 | < 1 GeV 2 . (102) 

This prediction is based on the EVMP HERA data, but the uncertainty in normalization 
on "standard candles" can be larger and reach a factor like *jj [85] . 

Then we can calculate also di-jet and di-gamma production. Results are presented in 
Figs. [121131 For the case of di-gamma production we can assume higher rate (as in CDF) 
due to nonperturbative effects. The latest LHC results |JT] give the upper bound for the 
EDDE di-gamma production at y/s = 7 TeV 

Op+p^p* +77+P* < 1-18 pb, E Til > 5.5 GeV, |r/ 7 | < 2.5, (103) 

and no particles in the region |r/ 7 | < 5.2. 

Now there are preliminary results from LHCb collaboration on exclusive Xc,o and J/ty 
production at ^/slhc = 7 TeV [12] : 




(104) 
(105) 



GeV 



Figure 12: EDDE di-jet production rate at LHC versus lower cut on the transverse energy 
of a jet for y/s = 8 TeV, \r]j et \ < 2.5 plus cuts f)102p . 
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Figure 13: EDDE di-gamma production rate at LHC versus lower cut on the transverse 
energy of a photon for = 8 TeV, |t/ 7 | < 2.5 plus cuts 



Our predictions are 



da 



dy 

excl, th 

Xc,0 

excl, th 
J/* 



dy 



a 



y=o 

excl, th 
J/* 



20 ± 5 nb, 
212 ± 53 nb, 

Clhc ■ (Tj+p-tj/y+piWo) = 7.06 ± 0.91 nb, 
76.3 ± 19.1 nb, 



(106) 
(107) 
(108) 
(109) 



where W = \Jmj/^^SLHC — 147 GeV and Clhc — 6-6 ■ 10 -5 . These results are close to 
the experimental data. 

Results for other interesting processes can be obtained by the use of master formu- 
lae (E]),®. 



4 Conclusions and discussions. 

In this article we consider the general framework for the exclusive double diffraction. In 
the present state of investigations we can describe it by the use of the prescription depicted 
in Fig. [5j 

• for high invariant masses (M ^> 1 GeV) we have to use a perturbative approach 
for the "bare" amplitude and calculate the gluon loop integral of the product of 
diffractive proton-gluon amplitudes T D , gluon-gluon fusion vertex and Sudakov-like 
suppression factor (Fig. |5^); 

• for intermediate invariant masses (M ~ 3 — > 10 GeV like in \ c production) we 
can apply the NRQCD approach to the QQ-quarkonium vertex and the gamma-like 
approximation for a Pomeron [83J; 
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• for low invariant masses (M ~ 1 GeV) we have only the general representation (189]) - 
( 19 2 p for the Pomeron-Pomeron fusion vertex. We can apply also some nonpertur- 
bative approaches like [54j-[23j. 

After the calculation of the "bare" amplitude we have to take into account different 
rescattering corrections ("soft survival probability"). 

Let us note that in this paper we use the old version of 3-Pomeron model only as an 
example. After the comparison with the latest TOTEM data [86], [87] it shows discrepancy 
between the prediction and data points (as other popular models [88]), and we have to 
some update our approach to better fit the data. 

The algorithm for EDDE calculations is far from an ideal one. This is because of our 
large gaps in understanding of the hadronic diffraction. There are different approaches 
at different scales. Also there are many models for diffractive amplitudes that are not 
always adequate (see [88] and references therein). This difficult task needs futher intensive 
investigations, and EDDE is the most powerful tool in this scope of activity. 
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